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Abstract 

In Split Supersymmetry scenarios the possibility of having a very heavy gravitino 
opens the door to alleviate or completely solve the worrisome "gravitino problem" 
in the context of supersymmetric baryogenesis models. Here we assume that the 
gravitino may indeed be heavy and that Majorana masses for neutrinos are forbidden 
as well as direct Higgs Yukawa couplings between left and right handed neutrinos. We 
investigate the viability of the mechansim known as Dirac leptogenesis (or neutrino- 
genesis), both in solving the baryogenesis puzzle and explaining the observed neutrino 
sector phenomenology. To successfully address these issues, the scenario requires the 
introduction of at least two new heavy fields. If a hierarchy among these new fields is 
introduced, and some reasonable stipulations are made on the couplings that appear in 
the superpotential, it becomes a generic feature to obtain the observed large lepton 
mixing angles. We show that in this case, it is possible simultaneously to obtain 
both the correct neutrino phenomenology and enough baryon number, making thermal 
Dirac neutrinogenesis viable. However, due to cosmological constraints, its ability to 
satisfy these constraints depends nontrivially on model parameters of the overall theory, 
particularly the gravitino mass. Split supersymmetry with 10^ GeV < »Ti3/2 ^ 10^° GeV 
emerges as a "natural habitat" for thermal Dirac neutrinogenesis. 



1 Introduction 



The phenomenology of split supersymmetry models ^ |2] in which a hierarchy is gen- 
erated between the gaugino masses and the masses of the scalar sparticles has received a 
great deal of attention in recent times. The main advantage of such scenarios is that they 
circumvent a wide variety of data pressures on theories with supersymmetry breaking on a 
lower scale that arise from potentially dangerous radiative corrections involving TeV-scale 
scalar sparticles for the Higgs Mass, flavor-violating effects, etc. This is done by making 
the scalar masses heavy, while keeping the gaugino masses at the TeV scale or below to 
constitute the dark matter and to preserve gauge coupling unification. 

In one particularly simple scenario, which we refer to as loop-split supersymmetry (after 
the PeV-scale supersymmetry of ^3_), supersymmetry is broken at an intermediate scale, 
around ~ 10^ — 10'' GeV, and all scalars in the theory, (with the exception of one light Higgs 
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particle) are given masses around the PeV scale while the gauginos acquire masses at the TeV 
scale [4, . The preferred method of achieving this hierarchy is by invoking anomaly mediation 
in the gaugino sector (but not the scalar sector) . This can be arranged by charging the chiral 
supermultiplet X responsible for the transmission of supersymmetry effects under some 
symmetry, to the effect that the term that normally provides the dominant contribution 
to the gaugino mass is forbidden by gauge invariance. The dominant contribution to the 
gaugino masses now arises only at the one-loop level |1] and is given by 



where the index A labels the gauge groups in the theory. Charging X will not affect the 
scalar masses, which are still manifestly gauge-invariant whether X is charged or not, and we 
obtain the desired hierarchy. The advantage of this scenario is its simplicity: no additional 
symmetries are required to keep the gaugino masses from being evolved up to the SUSY- 
breaking scale. The theory also imposes severe restrictions on both the identity (which 
dictates ought to be predominantly either Wino or Higgsino for strict one-loop AMSB) 
and mass of the lightest supersymmetric particle (LSP), and it has been shown that dark 
matter composed primarily of such particles could have observable consequences at the next 
generation of 7-ray telescopes [SUHIIZI- Also, in loop-split SUSY, characteristic signatures 
of gluino decays may be observable at colliders |H1 EI • 

As mentioned above, one of the reasons why split supersymmetry is such an attractive 
phenomenological model is that it provides a convenient way of circumventing data pressures 
on weak-scale superpartner masses. However, for a theory to describe the universe we live 
in, it is not enough that it evade all present experimental bounds: the theory must be 
cosmologically viable as well, in the sense that it does not disrupt big bang nucleosynthesis 
(BBN), that it is compatible with cosmic inflation, that it permits some mechanism by which 
baryogenesis may occur, etc. In this paper, we examine the cosmological ramifications of 
split supersymmetry, and in particular show that a viable model of baryogenesis can be 
realized via Dirac neutrinogenesis (i.e. Dirac leptogenesis) |lfl| lllj. This is not a trivial 
problem: it turns out that constraints arising from gravitino cosmology fix the reheating 
temperature associated with cosmic inflation to ~ 10^'^ GeV or lower for substantial regions 
of ~ iTT'LSP parameter space. This heightens the tensions already inherent among 
the model parameters (mass scales, couplings, etc.) of the theory. Data from neutrino 
oscillation experiments impose an additional battery of constraints any phenomenologically 
viable theory must satisfy. As a result, getting Dirac neutrinogenesis to work in split 
supersymmetry requires careful consideration of both cosmology and neutrino physics. 

Our aim in this paper is threefold. First, we examine the astrophysical and cosmological 
constraints on thermal leptogenesis models, especially those associated with baryogenesis 
and gravitino cosmology. Second, we investigate the neutrino spectrum constraints imposed 
on Dirac neutrinogenesis models by neutrino oscillation experiments. Third, we solve the 
full system of Boltzmann equations for Dirac neutrinogenesis and assess its viability — i.e. 
its ability to satisfy all aforementioned constraints. 




(1) 



2 



2 Dirac Neutrinogenesis 



The universe we live in is manifestly asymmetric between baryons and antibaryons. 
This statement can be quantified by introducing a parameter rj, defined as rj = nB/riy. 
Here ub = nb — n-^, where reft and re^ are the baryon density and antibaryon density of our 
universe, respectively, and is the present number density of photons. The value of r] has 
recently been measured with great precision by WMAP jl4j to be within the range 

r/ = (6.1 ±0.3) X 10-^°. (2) 

In the standard cosmology, which is symmetric with respect to baryons and antibaryons, 
one would expect to find = 0. As this is not to be the case, we must find a method of 
baryogenesis - the generic term for a process through which a baryon asymmetry might 
evolve in the early universe - to account for this. A set of generic criteria required for any 
successful baryogenesis scenario were first established by Sakharov ^21- first, there must be 
baryon number B violation; second, there must be C and CP violation; and third, there 
must be some departure from thermal equilibrium. If any one of these conditions is not 
met, baryogenesis fails^. 

A variety of viable baryogenesis models exist, including electroweak baryogenesis, in 
which the CP-violation occurs at a bubble wall, or phase boundary, and Affleck-Dine 
baryogenesis |18j . in which the baryon asymmetry is generated by moduli fields charged 
under B — L. In this paper, we will focus on models which achieve baryogenesis through 
a framework known as leptogenesis \W\ I2()j. In this scenario, decays of heavy particles in 
the early universe which violate both CP and lepton number L produce an initial lepton 
asymmetry, which is then converted to a nonzero baryon asymmetry by sphaleron processes 
associated with the SU(2) electroweak anomaly Leptogenesis is a particularly attractive 
model because in addition to its ability to yield a realistic value for rj [211, '^^^ ^is° explain 
why the standard model neutrinos have small but nonzero masses. In its most common form, 
which we will call Majorana leptogenesis, the heavy particle whose decays violate L is taken 
to be the right handed neutrino, which being a gauge singlet, may be given a large Majorana 
mass Mjy^. Since the most general realizable superpotential 

W B yLNHu + M^^NN (3) 

will also contain a term which gives rise to a Dirac mass mo = sin/3, where v is the 
standard model Higgs VEV, the neutrino mass matrix will contain small off-diagonal terms 
mixing vi and v^- When this matrix is diagonalized, the resulting mass spectrum contains 
three light neutrinos with masses 

rriu = mDj^m^jj, (4) 

which are identified with the standard model neutrinos, as well as three heavy neutrinos with 
masses ~ M^,^. This see-saw mechanism |21| is an example of one of the strongest assets 

^Actually the third condition is only required in theories where the Hamihonian preserves CPT. For 
a theory that is not CPT- invariant, such as the spontaneous baryogenesis of ITBl . departure from thermal 
equilibrium is not necessary for baryogenesis. In the scenarios we are considering, we will assume that CPT 
is a good symmetry of the Hamiltonian and that some departure from thermal equilibrium is required. 
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Figure 1: A schematic representation, after ^U], of the evolution of baryon number B 
(vertical axis) and lepton number L (horizontal axis) in Dirac and Majorana leptogenesis. 
In Dirac neutrinogenesis (left panel), the evolution of B and Ltot proceeds in three steps: 
first, two stores of lepton number (stored in left-handed neutrinos) and L^^ (stored 
in right-handed neutrinos) are produced during heavy particle decays; second, sphaleron 
processes (which act along lines of constant B — L) mix and B while leaving L^y^ alone; 
third, after sphaleron interactions have effectively shut off, equilibration between Li and 
Lu^ results in a net positive B and L for the universe. In Majorana leptogenesis (right-hand 
panel), only one store of lepton number is created, and the result is a universe with negative 
Ltot and positive B. 



of leptogenesis scenarios: they all provide some explanation for the lightness of neutrino 
masses, in addition to explaining the origin of the observed baryon asymmetry. 

Recently, another promising leptogenesis model, which is often called Dirac leptogenesis 
or Dirac neutrinogenesis jlUUllj . has received some attention. In this scenario, an additional 
symmetry is introduced, and charges are assigned under this new symmetry in a manner 
which forbids both the Majorana and Dirac mass terms appearing in (j^j) and additional 
heavy, vector-like fields are introduced whose decays will violate CP. These decays build 
up equal and opposite lepton asymmetries and Lj^^ in the left-handed lepton and right- 
handed neutrino sectors, while maintaining an overall lepton number for the universe Ltot = 
Li + Ljy^ = 0.^ In the fermion sector these stores do not equilibrate due to the smallness 
of the neutrino Dirac mass term, which only appears in the low-energy effective theory 
suppressed by powers of the mass scales associated with the heavy vector-like fields. The 
electroweak sphaleron processes which convert L^ into a baryon asymmetry B effectively 
shut off before L^ and Ly^ have a chance to equilibrate. The result, depicted in fig. is 
that the universe ends up with a net positive lepton number as well as baryon number. 

We will consider a model with the same field content as the one presented in In 
addition to the usual quark and lepton supermultiplets (of which the left-handed lepton 
multiplet L and the Higgs multiplets Hu and Hd will be pertinent to leptogenesis), we 
introduce a right-handed neutrino superfield A^, an exotic chiral multiplet and a number 

of vector-like pairs of chiral multiplets $i and $j (the precise value of is unspecified: 
any choice of > 2 is allowed from a baryogenesis standpoint). The charge assignments 
under the additional symmetry, whatever it may be, are to be arranged such that the most 

^Unlike in Majorana leptogenesis, no explicit lepton-number-violating terms are present here. 
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general superpotential that can be written is^ 

W 3 XiaNc^^iHu + hiaLa^iX + M$,$i¥i + fiHuHd, (5) 

where a is a family index, AjQ, and hia are Yukawa couplings, Mj are (supersymmetry- 
respecting) mass terms for the $i fields, and ^ is the usual Higgs mass parameter. The 
particular symmetry introduced and charge configuration employed in arriving at this 
superpotential is not of particular relevance to us, nor is it the aim of this paper to explore 
the model-building possibilities afforded by different such configurations, but an example 
(taken from of one that works is provided in tabled Mi, Xia, and hia may in general 
be complex. Both the scalar and fermionic components of the ^ and $ multiplets, which 
we denote by (j), (/>, ip^, and ip-^ will play the role that the i'r play in Majorana leptogenesis. 
As for neutrino masses, they will appear in the low-energy effective superpotential obtained 
by integrating out ^i and 

^eff B ^^xLHuN + fiHuHd. (6) 

If we arrange for the scalar component of x to acquire a VEV, then the first term in Weff 
translates into a neutrino Dirac mass 



This can be attained via an O'Raifeartaigh model of the type employed in in which the 
F-term of x acquires a large VEV (F) ~ 7713/2 -^p, and supergravity effects give rise to a 
nonzero VEV (x) — 167rr7i3/2«^~^ for the scalar component of X; where k is an undetermined 
dimensionless coupling constant. Requiring that (x) ^ where denotes the 

lightest of the M$. , and that the Xia and hia are 0(1) or smaller, this setup yields small yet 
nonzero neutrino masses without the aid of the see-saw mechanism. Furthermore, because 
of the simple structure of the neutrino mass matrix given in ((T)), Dirac neutrinogenesis can, 
under certain conditions, yield interesting predictions about the mass hierarchy among the 
standard model neutrinos. 

Just like in Majorana Leptogenesis, some intermediate scale is required to set the mass 
scale of the heavy fields responsible for the small neutrino masses and for the lepton number 
generation. In our case, we are less constrained since the neutrino masses are set not only 
by the masses M$. of the heavy fields, but also by the VEV (x) of the new field x- This 
leaves us more freedom to choose the mass scale of the heavy fields, but can also be thought 
of as a drawback since we have lost some predictiveness. One suggestive idea would be 
to make use of the x field and the additional symmetry to generate the supersymmetric fi 
term. The effective superpotential could then look like 



^eff 3 ^^^xLH^N + xHuHd. (8) 



^In engineering the charge assignments that lead to such a superpotential, one must of course worry 
about preserving gauge coupling unification, making sure the theory is free of gauge theory anomalies, etc. 
Since we are primarily interested in examining the phenomenology associated with such theories, we will 
not explicitly address these concerns, and appeal to high-scale physics (e.g. additional heavy fields, the 
Green-Schwarz mechanism |T2]) to resolve these issues. 
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Field U{1)l U{1)n SU{2) U{1)y 



N 


-1 


+1 


1 


L 


+1 





2 


Hu 








2 


Hd 








2 


<P 


+1 


-1 


2 


<t> 


-1 


+1 


2 


X 





-1 


1 



Table 1: One possible set of charge assignments, taken from that leads to the Dirac 
neutrinogenesis superpotential given in 0. Here the additional symmetry employed is a 
f7(l) (which may in principle be either global or local). Only the charges of the fields 
relevant to leptogenesis, which include the Higgs doublets Hu and Hd, the left-handed 
lepton superfield L, the right-handed neutrino superfield A^, the heavy fields $ and 
and the additional field x, have been included. Here, f7(l)L, SU{2), and U{1)y respectively 
denote lepton number, SU{2), and U{1) hypercharge quantum numbers, and U{l)]\f denotes 
the charge under the additional U{1). 

where is an 0{1) coupling constant, and where the Higgses should now be charged under 
the new hidden sector symmetry (and therefore enlarging the list of potential problems of 
the model, like cancelation of anomalies, etc.). In this situation, we have ^ = Y^{x) and 
any observation or limits on Higgsino dark matter would directly constrain the VEV {x)- 
(Split) supersymmetry might also give us some ideas as to how to relate the heavy fields 
to the SUSY breaking scale. We will come back to some of these issues when we deal in 
detail with neutrino phenomenology in section^ 

If Uia and hia contain nontrivial, CP-violating phases, a net CP asymmetry will be 
generated as the component fields of <I>i and <&i decay, the leading contribution to which 
results from the interference of tree-level and one-loop-level diagrams. Those relevant to (j) 
and (f) decay are shown in fig. The fermion fields and t/'^^ undergo similar decays, 
and stores of lepton number are built up in the lepton fields i^r and i, and in the associated 
slepton fields Dr, L The CP amplitudes (and resultant CP-asymmetries) generated from 
the decays of the charged fields in the fields charged under SU{2) will be the same as those 
generated from the decays of the neutral fields, since SU{2) is unbroken. Additionally, 
in the approximation of unbroken supersymmetry, the amplitudes (and resultant CP- 
asymmetries) in the fermion case will be the same as those for the scalar case. We will 
assume that any lepton number stored in the slepton sector rapidly equilibrates through 
the large (P^) term, and thus we need only pay attention to the lepton number stored in 
the lepton sector. 

We will begin by examining a toy model in which there are only two sets of ^ and 
$, the minimum number required for CP-violation to take place, and make the further 
simplifying assumption that the masses M$. are all much larger than /x. In this case |23) . 
one may parameterize the associated lepton number violation by defining a single decay 
asymmetry e, which represents the amount of lepton number generated in any particular 

*Here, we have included loop contributions involving the scalar component of Hd, which start to become 
important when ^ is of the same order as the Afs^ 
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Figure 2: Diagrams that give the leading contribution to the CP asymmetry from decays 
of the scalar fields (f)i and 0^. Similar CP asymmetries are generated during the decay of 
the fermionic fields and ifj— . 



lepton-number-carrying species by the decay of a single heavy particle. This implies the 
relations 
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where F^) is the total decay width of any of the heavy fields in the $1 or $1 supermulti- 
plets, and we have used the superfield notation for $1, N , etc. to indicate that the decay 
asymmetries are the same for all supersymmetrizations of the diagrams in fig. [21 Explicit 
calculation of F/j and e yields 



(13) 



and 



Im(At„A2a^^^%M^,M|J 
47r(|M$,|2-|M<,J2)(|Ai^|2 + |/,^^|2)' 



(14) 



where in both equations, a sum over the repeated indices a, (3 and 7 is assumed. Defining 
5 = |M$J/|Af$2|, e can be expressed in the more revealing form 



Im(At,A2^/i*^/i2/3e»^) / 8 
47r(|Ai^|2 + |/ii^|2) Vl-52 



(15) 



where if) is the relative phase between and M^^ . This tells us that for small values of 5^ 
the final baryon-to-photon ratio will be approximately proportional to 5. Since Ha and n/j^ 
have equal and opposite charges under the global U{1)l symmetry, the individual and 
Lyj^ lepton numbers respectively stored in left-handed leptons and right-handed neutrinos 
will likewise be equal and opposite. Thus no net lepton number is produced by the decays 
of ^> and the generation of both Ltot 7^ and B ^ occurs via electroweak sphaleron 
processes, which conserve B — but violate B + L^. As these are associated with the 
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SU{2) X U{1)y electroweak anomaly, they will only affect Li and not L^j^, and thus create 
both a nonzero value for B and a disparity between the two stores of lepton number (see 

fig. CD. 

In order for these two stores of lepton number not to be equilibrated away through the 
effective Higgs coupling, the equilibration rate must not become significant compared to the 
expansion rate of the universe until after the electroweak phase transition Tc, at which point 
sphaleron interactions have effectively turned off. This equilibration rate may be estimated 
on dimensional grounds to be 



^^1 172 9 T, (16) 



where g is an 0(1) gauge or top Yukawa coupling and T is temperature. The expansion rate 
may be expressed by the Hubble parameter H = l.Gdgl^'^T'^ /Mp, where is the number of 
interacting degrees of freedom (in the minimal supersymmetric model (MSSM) during the 
baryogenesis epoch, 205). Requiring that Tf,q < H ioi T > Tc leads to the condition 

\3m. < 10-8, (17) 

which can be translated into a limit on the neutrino mass using equation (f7|): 

< 1.74 sin /3 keV. (18) 

This condition has to be met in any case because of cosmological constraints on the absolute 
neutrino mass, and therefore left-right equilibration would happen well after sphalerons have 
shut off, and Dirac neutrinogenesis can then yield an appropriate value for r]. 

As we mentioned earlier, an important advantage of Dirac leptogenesis is its versatility: 
there is nothing in the theory that fixes any of the parameters M$^, M^j, Aja, hia-, and (x) 
to any physical scale. There are, however numerous constraints on these parameters. First, 
the neutrino masses must be appropriately small. Second, we must satisfy the Sakharov 
criterion that the abundances of (pi and (pi depart from their equilibrium values. This 
occurs when F/j is slower than the rate of the expansion of the universe at the temperature 
T ~ when (pi can no longer be treated as effectively massless and its abundance begins 
to fall off, or in other words 

9.97 • 10-^^/^:^ Ed^i-I' + < 1. (19) 



This constraint also favors small couplings and large M$^. However, in addition to these two 
requirements, we must ensure that the present value of r] satisfy the bounds in equation ((2)). 

To derive the precise limits this requirement places on the model parameters, one must 
solve the full system of Boltzmann equations, which we do in section|SJ For present purposes, 
a rough estimate can be made using the "drift and decay" approximation 24 , in which we 
assume that the (p decays occur well out of equilibrium and that the effects of inverse decays 
and 2 2 processes where AL^ 7^ are negligible. Including contributions from both scalar 
and fermion decays, this gives the result 

Le = ^ = -^K2{1) = 7.32 x 10'^, (20) 

s ■K'^g^ 
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where is the initial number density of and K2 (x) denotes the modified Bessel function 
of the second kind, evaluated at x. Since the final baryon-to-entropy ratio B (related to rj 
hy B = rj/7.04) generated by sphaleron processes will be on the same order (B ~ 0.35L^), 
this can serve as a rough estimate for B. Thus even if equation H19|) is satisfied, the final 
baryon-to-entropy ratio of the universe will be proportional to e; and from equation (|14j) . 
we see that for e to be large, either the couplings must be large or the splitting between 
and M<j,2 must be small. Thus there are tensions among the model parameters in Dirac 
leptogenesis, but they can be reconciled without too much difficulty. 

The real tensions among M,j,j, M^^, Aja, hi^, and (x) are not those inherent in the Dirac 
neutrinogenesis framework, however, but those that arise when we demand that the model 
respect the full battery of additional constraints from cosmology and from neutrino physics. 
We now turn to address these constraints and their implications for Dirac neutrinogenesis. 

3 Astrophysical Constraints 

Because the particle spectrum of models with high-scale supersymmetry breaking may 
differ significantly from that of models with weak-scale supersymmetry breaking, it is 
necessary to investigate whether any such alterations will affect the astrophysical constraints 
(from BEN, inflation, etc.) the theory must satisfy. In addition, Dirac neutrinogenesis 
introduces three additional light, sterile neutrino fields URa, which could affect BEN. We 
must make certain that neither of these considerations presents an insurmountable problem 
for our model. 

We will first address the issue of the additional neutrino species. It has been shown 
that the presence of the three additional light neutrino fields i^Ra at the time of EBN 
does not violate the bound [2^ on the number of additional light neutrino species Ni, < 
0.3, since their contribution to Ni, is suppressed by a large entropy factor. Still, there 
are other potential cosmological problems that might arise for split supersymmetry: in 
particular, constraints on the reheating temperature Tr associated with cosmic inflation 
become constraints on in thermal leptogenesis models. We must ensure that such 

constraints do not make Dirac neutrinogenesis unworkable. 

The most severe constraints on Tr are rooted in gravitino cosmology. There are two 
distinct varieties of gravitino problem that must be addressed: first, late gravitino decays can 
disrupt EEN by releasing energy in the form of photons and other energetic particles into the 
system; second, as the gravitino decays predominately into the LSP and its standard model 
superpartner, the potentially large, non-thermal population of stable particles produced in 
this manner could overdose the universe — or if the right amount is produced, could make 
up the majority of cold dark matter (CDM). Eoth of these issues are contingent on the 
gravitino lifetime T3/2, which for cases where J?> rnisp, is approximated by |22] 




(21) 



Careful analysis of the EEN constraints |27] reveals that unless ^ 10^ GeV, Tr cannot 
be greater than around 10® GeV, and for ^ 5 x 10^ GeV, cannot exceed 10® GeV. In 
models where m3/2 is at the PeV scale, however, r3/2 ~ 10~^ s, which implies that gravitinos 
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produced in the thermal bath decay long before the BBN epoch (at tuniverse ~ 1 s), and 
thus there is no gravitino problem of the former type for models with 77^3/2 at or above 
these scales. However, since a weakly interacting LSP decouples on a timescale i/ ~ 10~^^ 
s, it will have long frozen out by the time gravitino decay occurs unless is larger than 
around 10^ GeV; hence in theories with smaller gravitino masses (including simple PeV- 
scale supersymmetry with anomaly- mediated gaugino masses), LSPs produced by gravitino 
decay will be unable to thermalize and the latter type of gravitino issue cannot be ignored. 

In order to avoid any complications from late gravitino decay, we require not only that 
the LSP not overdose the universe, but that its surviving relic density ^lsp must be less 
than (or ideally, if the LSP is to constitute the majority of cold dark matter, equal to) the 
relic density of CDM as measured by WMAP |28] . 



OcDM/i^ = 0.11 ± 0.01 (WMAP 68% C.L.). (22) 

In general, ^lsp will have both a thermal and a non-thermal component, so that ^lsp = 
^LSP + ^LSP- "^^^ thermal component ^Igp may be ascertained by solving the relevant 
set of Boltzmann equations for the LSP abundance at freeze-out. The results, for the case 
where the LSP is essentially either a pure Wino or Higgsino, are [2j 



(\M2\\ 



nl^ph^ = 0.02 (^TjT^J for Wino LSP (23) 
nllph^ = 0.09 (^^) for Higgsino LSP. (24) 



Now we turn to evaluating ^^gp- We begin by addressing the regime in which there is 
no significant reduction in ^^Jp from LSP annihilations. Assuming for the moment that 
the dominant contribution to the non-thermal relic abundance comes from late gravitino 
decays and that all the LSPs produced from such decays survive until present day, ^igp is 
given by 



qNT 
^^LSP 



ruLSPrnTi^ 



T^"^ Peril 



3/2(^3/2)) 



(25) 



where pcrit is the critical density of the universe, Tq is the present temperature of the uni- 
verse, and 5^/2(^3/2) is the number of gravitinos per co- moving volume at the characteristic 
temperature T3/2 at which the gravitino decays, which is given by [29] 



^3/2(^3/2) 



0.856 X 10 



-n ( Tp \ 
VlOiOGeVy 



1 - 0.0232 In 



VlOiOGeVy 



Substituting this into equation yields 



^LSP^^ 



(26) 



(27) 



In fig. 121 we plot the contours corresponding to the WMAP upper and lower bounds 
from equation (|22j) on the total LSP relic abundance, as well as the simple overclosure 
bound Q.isph? = 1, taking into account both thermal and non-thermal contributions, as a 
function of misp and Tp. The region above and to the right of the WMAP upper bound is 
excluded. In order for the LSP to make up a significant fraction of the dark matter, ^lsp 
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Figure 3: Contours of ^^sph"^ corresponding to the upper and lower bounds from WMAP 
(|22|) . as well as a simple overclosure bound, as a function of the LSP mass misp and the 
reheating temperature Tr associated with cosmic inflation. The requirement that ^^lsp not 
overdose the universe severely constrains Tr, and hence the temperature scale of thermal 
leptogenesis, for a theory like PeV-scale loop-split supersymmetry, in which the LSP is 
particularly heavy. 
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should fall within the narrow strip between the upper and lower bounds on ^cdm given 
in (|22|) . but if some other particle makes up the majority of CDM, then the entire region 
below and to the left of the WMAP upper bound is phenomenologically allowed. In order 
to achieve an appropriate LSP relic density by thermal means alone, it can be seen that Tr 
must be quite low-around 10^ GeV; if the majority of CDM is generated non-thermally, Tr 
may be raised a bit, but is still constrained to be below ~ 5 x 10^'' GeV. As we shall see in 
sectional Tr ~ 10^ GeV turns out to be problematic for Dirac neutrinogenesis (in terms of 
the final baryon-to-photon ratio generated) , largely due to the out-of equilibrium condition 
in (dJ, this implies that if we want to raise Tr above 10^ GeV and still have the LSP relic 
density dominate i^cDM, the dark matter must be essentially non-thermal in origin. 

We now turn to address the regime where LSP annihilations do play a role in reducing 
^igp, and thus the upper bound on Tr may be raised. This effect becomes important when 
m3/2 ^ misp- When it is taken into account [Hll, the non-thermal LSP relic density is 
modified to 

ll^gp - mm ni^gp ,}^iLsp I, i28j 



where ^igp'^ is the relic density given in equation H27j) . and Q^gp"'^"^ is the relic density 



obtained by solving the full system of Boltzmann equations for the LSP. For a Wino LSP, 

.NT{a 
''LSP 



i i r cp IS given by 



^hsp - ^-41X10 (1 + ^^)3/2 UoOGeV J llOOTevJ 

(^LSP\\' [^^lfrn^\\^ (29) 



V "^3/2 J J \ 3 V "I3/2 J 
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where xw = mw /rriLSP] for a Higgsino LSP, which annihilates far less efficiently, Q.^gp^ 
will be even higher. 

In fig. |11 we show the relationship between rj^^*^""^ and misp for several values of 
m3/2- To be of any benefit in reducing ^igp, fl^gp"'"'"'^ itself must not exceed the WMAP 
upper bound, which is included for reference. From this plot it is evident that annihilations 
are only effective in reducing the LSP relic abundance below this bound when 771,3/2 is 
much larger than misP- This is a problem for loop-split supersymmetry with m3/2 around 
the PeV scale, where the ratio of misp to m3/2 is explicitly determined by the anomaly- 
mediation relation It is not, however, a problem for more general split SUSY scenarios, 
in which the splitting between rriLSP and can be much larger |29| I45j . Furthermore, 
when m3/2 is increased beyond around 10*, r3/2 becomes short enough that gravitino decay 
occurs before LSP freeze-out, and ^ijp drops to zero regardless of what the ratio of misP 
to rn,3/2 is. 

10^ f 1 1 1 1 1 1 1 1 : 
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Figure 4: The variation of Q^^^*^™^ (the LSP relic density function when LSP annihilations 
are taken into account) with LSP mass for several different values of rn-3/2- When 7713/2 = 

10^ GeV, only for rriLSP ^ 500 GeV will Q.^'^p"'^'^^ not generate an LSP relic density that 
exceeds the WMAP upper bound on the CDM relic density, and thus only for such an LSP 
mass can the reheating temperature Tr be increased beyond the naive upper bound from 
fig. 01 For larger gravitino masses 7773/2 > 10^ GeV, annihilations effectively reduce ^cdm 
and allow for T/j to be increased above the naive upper bound. 



Since thermal leptogenesis requires that a thermal population of and $1 be generated 
during reheating after inflation, we require that < Tr; hence an upper bound on Tr 
translates into an upper bound on . This means that there is substantial model tension 
when the Tr is constrained to be at or below ~ 5 x lO^*' GeV. This is not in any way 
peculiar to Dirac leptogenesis either: in a Majorana leptogenesis model, the constraints still 
apply with the lightest right-handed neutrino mass M^^ in place of -/Vf$^ . Here, the out-of- 
equilibrium condition ((T^ demands that Ai^ and hia be at most O(10~^) for M<j,j ~ 10^*^ 
GeV; on the other hand, H14|) and H20|) imply that unless there is a large hierarchy among 
the Yukawa couplings to different sets of <I> and ^ (so that Aiq <C A2q for all a), or unless 
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and M^^ are essentially degenerate, making the Xia and hia any smaller than 0(10"^) 
will yield insufficient baryon number — and we haven't yet taken into account the effects of 
2 2 processes and inverse decays. In short, leptogenesis becomes difficult when there 
are upper bounds on Tr, which occurs when 7713/2 ^ 10^ GeV and LSP annihilations are 
ineffective. 

There are two methods for increasing the baryon density in models where is tightly 
constrained by bounds on the reheating temperature, as indicated in equation ()14() . One 
possibility, to which we have already alluded, is somehow to arrange a large hierarchy among 
the Yukawa couplings to different sets of $ and $ (i.e. to adjust the matrixes A and h). Since 
the couplings to $2 only influence the physics via their CP- violating phases in e p4() and via 
their contribution to the mass matrix 0, increasing them will cause no phenomenological 
problems as long as a realistic neutrino mass spectrum is obtained. The second possibility 
is to make and M<^^ very close together in order to achieve a resonance condition in 
e (i.e. to adjust (5, the ratio of these two quantities). Since the perturbation theory we 
have used in calculating equation ((TH) is good as long as the separation of and M^^ 
is substantially greater than the value of the off-diagonal elements Mij in the $ — $ mass 
mixing matrix induced at the one-loop level: 

M^j ~ g*g' {Nh, + )hj , (30) 

where gt and g'j represent the appropriate A and h, summed over the fermion family index, 
and lij is a numerical factor on the order of l/167r^ from the loop integral. Thus for small 
A and h, 6 can be set very close to one and M,^^ and M^^j may be very nearly degenerate. 
This allows for the possibility of resonant leptogenesis, as has been done in the Majorana 
leptogenesis case, which can be invoked to generate a large baryon number in cases where 
the bounds on Tr are more severe. As we shall see in the next section, scenarios with small 
S can naturally produce a realistic neutrino spectrum; thus we will restrict our attention to 
the case where 6 is small. 

While the problems that can arise for small gravitino masses have now been thoroughly 
addressed, the caveats associated with extremely large should also be mentioned. As 
has been shown in [H^, split supersymmetry models with a large hierarchy between the 
gravitino and gaugino masses can suffer from phenomenological problems associated with 
the overproduction of gluinos. While a precise ceiling for 7713/2 must wait until gluino decay 
is better understood, it is known that this ceiling falls somewhere in the 7773/2 — 10^" — 10^^ 
GeV range. There are also caveats associated with the gravitino-producing decays of scalar 
sparticles in models where one or more scalars has a mass larger than 7773/2 j32j . 

The message here is that cosmological considerations can make thermal Dirac neutrino- 
genesis substantially more difficult when 7773/2 < 10® GeV, as substantial tensions arise 
among the out-of-equilibrium decay criterion, overclosure bounds related to the reheating 
temperature Tr, the equation that determines the decay asymmetry e, etc. For light 
gravitinos (777,3/2 ^ 10^ GeV), BBN limits severely constrain Tr, and hence M$^. For slightly 
heavier gravitinos (10^ GeV < 7773/2 ^ 10® GeV), there are still constraints on Tr from 
nonthermal decays which are only alleviated (via LSP annihilations) when raisp "13/2- 
Since this is not the case in loop-split supersymmetry, where the ratio of rriLSP to 7773/2 is 
specified, making Dirac neutrinogenesis succeed is quite difficult in this framework. Still, if 
we allow for non-thermal generation of ^CDM (or for the possibility that the LSP is not a 
significant contributor to the dark matter density) , Dirac neutrinogenesis may still be viable 
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in loop-split SUSY (while there are apparent tensions in this case, they an be resolved by 
invoking some mechanism, such as resonant leptogenesis, for amplifying the lepton number 
produced). In more general split supersymmetry scenarios, however, Tr and may be 
increased to lO^'^ GeV or more, leptogenesis should proceed without too much trouble, and 
the possibility for thermal dark matter still exists. Whatever the case, it will be necessary 
to solve the full system of Boltzmann equations to ascertain whether or not there is any 
region of parameter space in which Dirac neutrinogenesis can be successful, which we do in 
section |SJ 



4 Neutrino Physics 



In addition to respecting constraints arising from cosmological considerations, in order 
to be phenomenologically viable, a given Dirac neutrinogenesis model must yield a neutrino 
spectrum that accords with current experimental constraints. The most stringent such 
constraints come from solar and atmospheric neutrino oscillation experiments, and place 
limits both on the mass splittings 

Ami,, = ml^ - ml^, (31) 

where the indices a and b label the different neutrino mass eigenstates, and on the mixing 
angles 6ab between these eigenstates. The primary connection between the latter and 
observable physics occurs through the leptonic mixing matrix 

Umns = U^''^U^'\ (32) 

where U^'^^ is the neutrino mixing matrix and U^^^ is the charged lepton mixing matrix. In 
the basis where the charged lepton mass matrix is diagonal, Umns is the neutrino mixing 
matrix and may be expressed in terms of the neutrino mixing angles 6ab as 

C12C13 S12C13 Si3e*'^C'P \ 

Umns = \ -S12C23 - ci2S23Si3e*'^^'' C12C23 - si2S23Si3e"*'^^'' S23C13 , (33) 

S12S23 - Cl2C23Sl3e"*'^C'P -C12S23 - Si2C23Sl3e~'^^P C23C13 / 

where Cab = cos 6ab, Sab = sin^ab) and 5cp is a C-P- violating phase. The present limits^ on 
the Am^jj and 6ab are [25 

sin^ (9i2 = 0.30;[};g^, sin^ 6*23 = 0.50l:[J:J^, sin^ ^13 < 0.031, 
Ami, = X 10-'eV2, \Aml^\ = (2.21°:^) x IQ-^eV^. (34) 

The smaller of the two mass splittings, Am^i, is to be identified with the Attiq obtained 
from solar neutrino data (the MSW-LMA solution); the larger, Am|,, with the Am^ from 
atmospheric neutrino data. The constraints in (|34|) imply a Umns niatrix with bounds 

- .86 .49 - 

\Umns\=\ .30 -.58 .40 -.68 .61 - .80 |, (35) 

- .46 .50 - 





^We do not take into account the LSND result which would require an extra neutrino mass eigenstate. 
In the case that forthcoming data from experiments such as MiniBooNE corroborate the LSND signal, it 
will be necessary to extend the neutrino content of our model. 
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but they leave the sign of the largest mass squared difference Am|^ undefined. This implies 
that the physical neutrino ^3 can be either the heaviest of the three mass eigenstates, i.e. 
mi < 771-2 "T-s (normal hierarchy) or the lightest, i.e. 7713 <C mi < m2 (inverted hierarchy). 

The unitary matrix Umns-, in the charged lepton basis, is responsible for diagonalizing 
the squared neutrino mass matrix , i.e. 

[mlY"^^ = U'^MNS rnvml Umns- 

We can therefore estimate the generic form of the neutrino mass matrix squared, since 
it has to be diagonalized by the Umns rnatrix. 

We should have for a Normal Hierarchy 

(m2)^-- ~ Ami I • 1 1 h (36) 




and for an Inverted Hierarchy 

{mlY^^^Amlil .111, (37) 




where the are small entries (not necessarily equal) compared to the 0(1) entries. In both 
cases, the order of magnitude of ^ is constrained by the ratio P32 = ^^31/^1^21, which, 
according to (|34j) . must respect the bounds 

20.0 < P23 < 39.7 . (38) 

We are interested in finding a simple connection between the heavy ^ fields and the 
neutrino sector, which does generically reproduce the observed neutrino spectrum and is 
also compatible with successful baryon number generation. 

Let us first consider the spectrum of a model with only one set of (j) and (f), or alter- 
natively, a theory in which <C M$- for all i > 1. In such a case, the mass matrix is 
proportional to an outer product of the family-space vectors \ia and thus its eigenvalues 
are 



m 



77i^2 = '^i^s = ^^-^la^la- (39) 



For each additional set of (f) and with a mass similar to M^^, an additional neutrino 
acquires a nonzero mass. Thus in "short-suited" models in which there are effectively only 
two sets of (j) and (p involved in determining the neutrino spectrum, one neutrino mass 
eigenstate is massless, and the other two neutrino masses are given by m^^ ^ ^^ii ^"^^ 
~ Am|]^. Furthermore, leptogenesis in such models is well-approximated by the toy 
model considered in section [2 in which there were only two sets of $ and 

From ((TJ, the neutrino mass-squared matrix is given by 



|2 

lo/3 



v{x)sm. f3\ ^ 



2 (or 3) 3 

E 

«J=1 7=1 



X*^Xjyh*^hji3 



(40) 
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As long as A and h are completely generic and there are at least three sets of (j) and <j), it is 
apparent that a matrix of this form can yield an arbitrary neutrino mass spectrum. On the 
one hand this is good, for it means that the theory is perfectly viable; on the other hand, 
this arbitrariness comes at the price of introducing many additional free parameters, whose 
relative values must be determined by some additional underlying physics. 

Guided by the idea of linking the neutrino data to the ^ fields sector, we find two 
simple conditions that generically reproduce the form of the normal hierarchy neutrino 
mass squared matrix 



Normal Hierarchical ansatz A 

A quick inspection of (|4U() shows that if the spectrum of heavy fields is hierarchical, 
i.e. if <C M$2(<C M.^.^)^ then we can very simply write the neutrino mass squared 
matrix to zero order m. 5 = M$j/M$2 as 



/ \2 Ai / '^^^'^ ^11^12 h\^hu . 

I"i.|' = (^(x)sm/3) ^ . |/ii2p hl^h^^ \+0{5), (41) 



where ^^i = Yj-^\^Ii? ■ 

Another hierarchical assumption, namely the requirement that the off-diagonal ele- 
ments of the coupling matrix h be much greater than the diagonal elements {hi2 ~ 
/ii3 = h ^ /ill), allows us to rewrite the previous equation as 



\m,\^ = {v{x)smpf^-^ 1-11 \ +Oi5), (42) 



£2 e 



where e = hu/h. 

This reproduces the structure required for a normal hierarchy neutrino sector. It is 
interesting that the two simple hierarchical requirements, ^ M<j>2 3 and hn <^ 
hi2 ~ /ii3 (with /ill 7^ 0), generically give rise to the correct neutrino phenomenology. 

Normal Hierarchical ansatz B 

The previous ansatz suggests that we look for some setup that would ensure that the 
diagonal element hn is small compared to the off diagonal terms hi2 and /ii3. An 
antisymmetric condition on the matrix h can do the job, but from eq. (|41() it becomes 
clear that we need to include a higher order term in the small parameter 5. 

When we assume that all the diagonal elements of h are zero (or extremely small), 
and the rest of elements in /i or A are of order 1, we obtain the structure 



(43) 



Again we find the generic structure required for a normal hierarchy, but now the value 
of 6 fixes the ratio y023- 




^Of course Mif,^ cannot be too large since then, as remarked earlier, we effectively generate the mass of 
just one physical neutrino leaving the other two massless. 
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Clearly, either ansatz A or B (or a combination of both) can explain the neutrino data 
and large mixing angles. Ideally, we would like to have a theoretical framework that would 
provide both the form of the superpotential required for Dirac neutrinogenesis © and the 
necessary flavor structure, like the ones presented in the previous examples. A careful 
construction of such a model is out of the scope of this paper, but we nevertheless would 
like to point at least in one interesting direction in which simple assumptions enable us to 
reproduce the needed conditions, while at the same time reducing the number of new free 
parameters. 

In many grand unified theories and models with non-Abelian flavor symmetries, the 
Yukawa matrices (and in general the flavor interactions) can be symmetric, antisymmetric, 
or both (see for example |35[ I36[ l37] and references therein). Operating in this paradigm, let 
us assume that the SM left-handed leptons have the same flavor charge as the heavy fields 

and the SM right-handed charged leptons have same flavor charge as the fields Upon 
breaking of the flavor symmetry, let us assume that the charged lepton Yukawa matrix is 
symmetric due to a symmetric flavon VEV configuration {Sa/s)- The corresponding effective 
superpotential is 

WDYi ^^HaL^ep + (5,^)l>,$^+h.c. (44) 

where Yi and Y$ are dimensionless couplings, and a and fi are flavor indices. The charged 
lepton Yukawa matrix {yi)a/3 and the mass matrix (M<j>)a^ of the heavy flelds ^ would be 
both symmetric and proportional 

{M^)ap = iy,)ap (45) 

This speciflc structure predicts exactly the mass spectrum for the (p flelds in terms of the 
flavor scale Mf, which is at the origin of the intermediate scale required for successful 
thermal Leptogenesis. 

Mf Mf Mf 
M<j, = rrie , -/Vf$„ = rrin and M^^ = nir (46) 

V V V 

If the flavor scale is of the order of some GUT scale Mf ~ 10^^ GeV, we would then expect 
M<j,j ~ 10^^ GeV, with 6 = M^^/M^^ = •nif.jvn^ ~ 5 x 10~^ being a small parameter. If 
one worries about the dependance on the Flavor symmetry scale it is also possible to link 
the intermediate scale to the large SUSY breaking scale MsusY- We can imagine that if 
all flelds in the theory are charged, the spurion fleld X of supersymmetry breaking might 
as well be charged, and perhaps it is charged under the same symmetry U{1)m responsible 
for obtaining the superpotential ©. Before any kind of breaking we can imagine the 
superpotential as 

W:^^{Yi HdL^ep + Y^X ^^^p) +h.c. (47) 

IVl F 

where again Yi and y$ are dimensionless couplings. Supersymmetry breaking effects can 

1 /2 

provide the fleld X with a VEV {Ax) ~ {m^/2Mpi) > and upon flavor symmetry breaking 
we could get the effective superpotential 



Weff D (yJa/3 [HdLaep + ^m^/2Mpi ^a^pj +h.c. (48) 

where we have assumed that the original constants Yi ~ 1$. In this situation, for example 
with m3/2 ~ 10^° GeV, the mass of the hghtest ^> fleld would be ~ 10^^ GeV. 
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Now, let us also assume that the flavon VEV {Saf3) is symmetric and that the coupling 
^Lx becomes antisymmetric upon flavor breaking, i.e. the superpotential can be written 
as 

W D {y,tyH,L^efs + {M^yl^^a^p + X^pN^^pH^ + /^^^^^^^'"i^L^x (49) 

where we have A^^ = (^^) and hf^^^ = (^^^), with the flavon acquiring an anti- 
symmetric VEV configuration and with the VEV configuration of being arbitrary in 
flavor space (it is not important to our purposes its flavor structure although probably in 
a specific model of flavor it would end up being some linear combination of symmetric and 
antisymmetric VEV configurations). Let us see what the implications of these ingredients 
are for our model. 

In general, the charged lepton Yukawa matrix (yja/s can be diagonalized by a biunitary 
transformation of the type 

[^ty^y=y(-i.a,)^ (50) 

but when yi is symmetric, this biunitary transformation takes the simpler form 

U^y,U=yf'^^l (51) 

If (M$) oc (yj, as in the setup described above (see equation (HHl)), then (M$) will be 
diagonalized by the same transformation. When the mass matrices for the charged leptons 
and the (p — (p system are simultaneously diagonalized (that is, when we go to the charged 
lepton basis), A and h transform as 

A' = U^XU h' = U^h-^^'U. (52) 

Transformations of this type preserve the antisymmetry of h (and A if also antisymmetric); 
thus in the charged lepton basis the matrix h remains antisymmetric and (M$) is real and 
diagonal; this reproduces the Hierarchical ansatz B described earlier, in which the diagonal 
elements of h are zero and it is the smallness of 5 the responsible for the large mixing angles 
observed in the lepton mixing matrix Umns- 

In the setup discussed above, (M$) and yi could be simultaneously diagonalized because 
(M$) oc yi, which fixes 6 = M<j,j/M$2 = rrie/m^. Such proportionality is not requires, 
however: any matrix of the form 

(M^) = Ayi + S/3X3, (53) 

where A and B are arbitrary constants and /sxs it the 3x3 identity matrix, can be 
diagonalized along with yi. Thus when B ^ 0, 6 can in principle take any value (as long as 
that value is consistent with the observational bounds on the neutrino spectrum). 

• Constrained Hierarchical Dirac Leptogenesis 

Motivated by the preceding remarks, we define Constrained Hierarchical Dirac Lep- 
togenesis (CHDL) as the phenomenological setup in which (M$) is real and diagonal, 
A and h are antisymmetric, and the smallness of 5 = M^^/M^^ explains the large 
mixing angles in the neutrino sector. We may parameterize A and h as 
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/ 1 as \ 
A = / -1 as 

\ -02 -as / 



/ 6i 62 \ 
h = f[-bi 63 , (54) 
V -b2 -h / 



which is convenient when the a and h are all roughly 0{\). Since the assumption of a 
hierarchy among the M$ . leads to the neutrino mass-squared matrix of the form (|41j) , 
we expect that 03 and 63, the effects of which show up only at the 0{5) level, will be 
less tightly constrained than the rest of the Cj and 6i, which contribute to the leading 
term. This is in fact the case: if ai, as, 61, or 62 deviates significantly from one, the 
neutrino spectrum cannot satisfy the constraints in (|34|) . It is therefore appropriate, 
since the value of / is unimportant as far as this set of constrains are concerned 
(any rescaling of / can be compensated for by a similar rescaling of {x)), to analyze 
constrained hierarchical models as functions of 03 and 63 alone. 



In fig. 121 we show the region of viability in — 63 space for two different values of 5: 
in the left-hand panel, we set 5 = me/m^ (and M^^/M^^ = m^/rrir) as required in the 
minimal version of CHDL discussed above; in the right-hand panel, we set 5 = 10~^. We 
consider a given combination of as and 63 to be phenomenologically viable if there is any 
combination of the remaining a^ and hi for which the combination simultaneously obeys all 
the neutrino oscillation constraints in (|34|) . This plot demonstrates two important features 
of the Yukawa matrices in CHDL: first, it is indeed possible to satisfy the neutrino oscillation 
constraints for 5 = m(,/m^] second, while ^3, like most of the other aj and ftj, is constrained 
to lie fairly close to 1, as can be quite large when 6 is small. 

In the same figure, we also show contours for the value of sin^is, the value of which 
will be measured or constrained in future neutrino experiments. It is seen that the value of 
sin ^13 increases with increased as until reaching its maximum experimental bound. 

Let us now take a moment to address how these results affect leptogenesis. Since we 
are assuming that » , M^^ , the formula ((T3)) for the decay asymmetry e tells us 
that lm.{X\^\2ah\ph2pMi^^M^^) will vanish (when it involves diagonal elements of A or h) 
unless (3 = a = 2>. This means that 

e (X Im (asas^s^s) -^—^ (55) 

in CHDL; the two panels in fig. [5] then show that for a given 5, the largest amount of left- 
handed lepton number (and therefore the largest baryon asymmetry) will be obtained 
when as = a'^°'^{5), where a^°'^(5) is the maximum possible value of as for a given 6 
consistent with neutrino masses and mixings. It is interesting to note that since the 
maximum experimental value of sin^is sets the value of a™"^^, then the maximum baryon 
asymmetry will be obtained when sin ^13 acquires its maximal experimental value, that we 
take to be sin^^^ ^13 = 0.031. Because S = me/m^ ~ 4.83 x 10~^ is fixed to be quite 
small and since e is approximately proportional to 6 for small 6, one might worry that 
such a small 5 might doom baryogenesis. However, as indicated in the left panel of fig. El 
a^"''^ {me/rrifj^) ~ 95, so that a hierarchy between couplings to different sets of ^ and $ is 
permitted, and the result is that e is only suppressed by an 0{1) numerical factor. The 
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Figure 5: Here the regions of 03 — 63 space (see equation for a description of the Yukawa- 
matrix parametrization used) for which all constraints on neutrino masses and mixings H34() 
are simultaneously satisfied for some combination of the remaining and bi are shown for 
two different values of 6. Additionally, contours depicting the ranges for S13 = sin ^13 (which 
depends primarily on 03, but varies slightly with the remaining and bi) are shown. In 
the left-hand panel, S = M^-^/M^^ = nie/m^; in the right-hand panel, 6 = 10~^. The plots 
reveal that while 63 is constrained, along with most of the other Ui and bi, to lie reasonably 
near 1, a^°'^ is on the order of 10^, which, when 03 ~ a™"^, results in an increased e and — at 
least in the drift-and-decay limit — the baryon-to-photon ratio r], making it easier to achieve 
a realistic value for r]. For smaller 6, a™'^^(5) is much lower, as will be rj. In each panel, the 
configuration that yields the greatest decay asymmetry e is marked with an arow. 

right panel of fig. [S] shows that a'^'^^{6) drops sharply as 6 increases and 0{6) corrections 
become important in m^, since 0^*^^(1/10) ~ 4.5. 

We are also assuming that we have maximal CP violation in the decays of the fields 
$ and $ (i.e. the overall phase in the product 020362^3 must be tt/2). In this case one 
can obtain a value for the effective CP violation in the lepton sector. This can be defined 
in a phase invariant way in terms of the quantity J = Im(C/i2?7|2t^23f^i3) EHl where 
Uij = (UMNs)ij- Taking for example the two points that will generate the most baryon 
number in the left and right handed panels of fig.[21(marked with an arrow) we find J ~ 0.034 
for the point in the left panel and J ~ 0.030 for the point in the right panel (since the 
maximal value for J is |C/i2f^22 ^23^^13 1) which for these points yields a value of J ~ .038, 
the CP- violation here is close to maximal). It could be interesting to do a more detailed 
study of the issue of linking more generally the effective CP violation in the lepton sector 
to the CP violation in the interactions of the heavy fields <I> and Nevertheless, we will 
not pursue this issue further in the present work. 

To summarize this section, we have shown that it is possible to construct a variety of 
Dirac neutrinogenesis scenarios (including the simple, theoretically motivated CHDL) that 
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are capable of satisfying current constraints on the light neutrino spectrum. In CHDL, which 
wiU become our primary focus from this point forward, the set of parameters relevant to 
neutrino physics and to leptogenesis is quite small, given the antisymmetry condition on the 
two coupling matrices A and h. The other important parameter is the ratio of masses 6 = 
M$^/M$2. It remains to be seen whether CHDL — or indeed any Dirac neutrinogenesis — 
model can also simultaneously satisfy the constraints from gravitino physics and yield a 
realistic baryon number for the universe. As we shall see in the next section, thermal Dirac 
neutrinogenesis will indeed turn out to be workable in a variety of split-supersymmetry 
models. 



5 Boltzmann Equations 



Let us now turn to the numerical calculation of r] and the solution of the full Boltzmann 
equations. The heavy fields aside, in Dirac neutrinogenesis there are six particle species 
charged under lepton number (z^/j, Dr, i, i, and the right-handed charged lepton and slepton 
fields cr and br), and thus six individual stores of lepton number to keep track of: Li, L,^^, 
L^, Lyj^, Lej^, and Lg^. However, since £, £, eji, and eji participate in SU{2) and/or U{1)y 
gauge interactions, which we assume to be sufficiently rapid (compared to other processes 
relevant to leptogenesis) that these species will always be in chemical equilibrium with one 
another, any lepton number stored in any one of them will be rapidly distributed among 
L£, L^, Le^, and Lg^ in proportion to the relative number of degrees of freedom of each 
respective field. Furthermore, since we are assuming that d. and vr equilibrate rapidly 
through a large {F^) term, ur will also be in equilibrium with these fields. Thus vr^ £, 
£, e/j, and cr can be viewed as one lepton sector with an aggregate lepton number Lagg- 
This leaves only the right-handed lepton field vr^ which does not participate in any of the 
aforementioned rapid interactions and thus retains its own distinct lepton number L^^. 

In the limit of rapid equilibration within the Lagg sector, the set of Boltzmann equations 
governing the evolution of lepton and baryon number simplifies considerably. In addition 
to the relations that describe the dynamics of the heavy fields, only three equations are 
required: one for Lagg, one for L^,^, and one for B. However, before we state these equations 
(a derivation of which is provided in Appendix^, it will be useful to digress for a moment 
and discuss the rates for the various processes involved. We then return to a discussion 
of the Boltzmann equations themselves, followed by a presentation of the results from our 
numerical computations. 

While there are a large number of 2 2 processes which transfer lepton number between 
different particle species, we will only need to calculate the rates for those which transfer 
it between Lagg and L^^. As for the rest, which collectively serve to equilibrate lepton 
number among the fields in the Lagg sector, it only matters that they are rapid compared to 
other rates in the problem. The diagrams for the s-channel 2 <-> 2 processes which shuffle 
lepton number between Lagg and L^^^ are pictured in fig. IHl In addition to these, there are 
contributions from the t-channel transforms (two per diagram) of these diagrams. In order 
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to evaluate diagrams containing virtual heavy fermions, we define a Dirac spinor 



where (V'3>i)a and are the Weyl spinor components of the and <I>j superfields. 

Numerical calculation of the thermally averaged cross-sections for the diagram pictured on 
the left in the top row of fig. IHl which involves two Yukawa-type couplings of a scalar to 
two fermions, yields {av)iap — cr) | Aja | | Aj^j | , where 

For the diagram on the right in the first row of fig. EI which include one Yukawa-type 
coupling and one trilinear scalar coupling proportional to M$-, the result is very nearly 

temperature independent and well approximated by (crv)ia/3 — fj-^^^^^IAiaPl/ij/jp, where 

<,1'™'-0.5X ' , (58) 

^^^^^ 

The two diagrams in the second row of fig. IHl which involve two Yukawa couplings and a 
mass insertion from the heavy fermions, yield this same contribution. These interactions 
dominate among 2 <-> 2 processes. The diagram in the bottom row of fig. El which involves 
a trilinear scalar coupling to the down-type Higgs, may be approximated by {crv)iaf3 — 

crl^'^^\Xia\'^\\i(i\'^, where 

where the constant of proportionality is 0(1), and is thus suppressed relative to the rate 
given in ()58(1 by ///M$-. Here, we will assume that // is several orders of magnitude smaller 
than all the M$. , and therefore the effect of these processes can be neglected. 

Taking into account contributions involving virtual fields in the $2 and $2 supermulti- 
plets, as well as those in <I>i and $1, the total interconversion rate between Lagg and Ly^^ 
is 



r2^2 ^ 3n^J^J^((Tf'')|Ai,|2|Ai^|2 + a?^)|A2a|'|A2/3p 

a 13 

+ 9n,5;j;(ar^^)|Ai„n/.i,P + ar''^|A2«|>2,l^) 

a (3 

+ 3n,^j;(a(^'^)|Ai„|2|A,^|2 + af'^)|A2„P|A2,p) (60) 

a f3 

where we have assumed an equilibrium number density for all non-leptonic light species 
(e.g. Hu, x) involved. 
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Figure 6: Diagrams for 2 <-> 2 s-channel processes which transfer lepton number between 
Lagg (the aggregate lepton number in the sector comprising the fields £, £, Dr, cr, and en, 
which are assumed to be in chemical equilibrium with one another due to rapid gauge and 
{Fy^-texm. equilibration interactions); and L^^ (the lepton number stored in right handed 
neutrinos). The two t-channel interactions associated with each diagram appearing above 
must also be included in calculating the full thermally averaged cross-section. 

In addition to the 2 2 processes acting, we must also discuss the rate associated with 
sphaleron processes, which represent saddle-point transitions between distinct electroweak 
vacua possessing different values of B and Lagg, for these will be responsible for the 
conversion of lepton number to baryon number. As these processes are mediated by the 
SU{2) electroweak anomaly, they only affect Lg (and hence Lagg) and not L^,^. We will 
primarily be concerned with the high temperature sphaleron interaction rate, by which we 
mean the interaction rate at temperatures T ^ Tc, where Tc is the temperature at the weak 
scale. This rate takes to form 

'^Sph = cTal, (61) 

where a2 = ^l/^'"" ~ 1/30 and c is an 0(1) constant, whose value calculations j42l 1431 HI] 
place very close to 1. 

We are now ready to write the full set of Boltzmann equations governing the evolution 
of baryon number B in the early universe. We will need to include include equations for 
Lagg, L^y^, B, as Well as ones for the lepton numbers stored in the heavy fields (pi and (pi, 
which we respectively dub and L^— , and the number-density-to-entropy ratio of one of 
these heavy fields or their conjugates (we choose (p^), which we call Y^^- In terms of the 
variable z = M^^/T, they are 
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dL 



dB 

dz 



z 



dz 



z 



{TR)L.-2Lagg[{TD),o+{TL) 



(62) 



+ {Luji - -Lagg){T2^2) 



dL 



dz 



-;;Lagg){^ 2^2) 



dY^ 



dz 
dL^ 

dz 
dLA,_ 



z 
z 



-{Td){Y^^ - Y;1) + ]^Lagg{TL),u + (rii) 

-{Td)L^^ + 2Lagg{TL)iD + 2Ly^{TR)j 



-{^d)L^_ + 2Lagg{TD) 



dz H{M^^] 
Here, the inverse decay rates (r/))^^, (ri^)^^, and {Tji)jjy are defined by 



(63) 

(64) 
(65) 
(66) 
(67) 



eq 


^K,{z)^ 


7 rij ' 




eg 


(K,{z) 


7 n-y 


\K2{z) 


eg 

_ 0* 


(Ki{z) 




\K2{z) 



r^,, and 



(69) 
(70) 



where F^j is the total decay width of 4>i given in equation ((T^ . n^*^ is the equilibrium 
number density of (jii (and (/>!, etc.), and the quantities Tl and Tr represent the partial 
decay widths for vrH (or (j) — > vrH) and cj) ^ Ix (or (/> ix)- The ratio of modified 
Bessel functions Ki{z) and K2{z) appearing in these rates is a result of averaging over time- 
dilation factors (see ()85p in Appendix^. Explicit definitions for these quantities, along 
with a derivation of the Boltzmann equations themselves, are provided in Appendix 1X1 

We solve this system of equations numerically as a function of the input parameters 6 
and the reparameterized coupling strength /, defined in eq. (|54|) . and present our results in 
fig. [71 Here, the regions of / — (5 parameter space in which the final value of r] generated falls 
within the WMAP-allowed range given in (j2j, appear as thin 'ribbons', each corresponding 
to a different value of . In fig. [71 the effects of the processes detailed above are apparent, 
and certainly nontrivial. Physically, these effects can be interpreted as follows: increasing 
the strength of the neutrino-sector couplings (here parameterized by /) increases Td, which 
in turn increases the initial value of B; however, from equation (|6U() . increasing / also 
increases the rates for the 2 <-> 2 processes which shuffie lepton number back and forth 
between Lagg and Lr. Furthermore, it increases the rate for inverse decays. This allows 
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two possibilities for generating a realistic final value for B. In the first case, where / is small, 
the initial baryon number produced by 4> and (j) decays is approximately within the range 
allowed by WMAP, and 2 <-> 2 and inverse decay processes are so slow as to be negligible; 
this is the "drift-and-decay limit" of equation (|2U() . In the second case, where / is large, a 
surfeit of baryon number will initially be produced, but these processes, which occur more 
rapidly for larger /, subsequently reduce i? to a phenomenologically acceptable level; this 
we refer to as the " strong- washout regime". These two possibilities are shown in the two 
panels of fig. |H1 in which the dynamical evolution of B, Lagg and other relevant quantities 
has been plotted, for 5 = 4.83 x 10^^ and M$ = 10^^ GeV. In fig. [3 the two regimes are 
represented respectively by the lower and upper portions of each ribbon. 

Having discussed the general effects of the washout processes described above as a group, 
it is also important to address their characteristics relative to one another. Inverse decays 
dominate over 2 <-> 2 processes only for a brief period, where 1 ^ ^ ^ 50, but during this 
period they are extremely effective in reducing lepton number, and in fact are the primary 
factor in determining the final value of ??. For larger z, until they freeze out, the 2 <-> 2 
interactions dominate and further reduce Lagg and L^^ (and consequently B). It should be 
noted, however, that the total B — Ltot number of the universe is manifestly conserved by 
the Boltzmann equations (|62|) - (|67j) (the sum of the rates for the various lepton numbers 
involved is zero), and since we began with B — Ltot = 0, we end up with B — Ltot = in 
any case. 

For a given and 5, then, there are at most two points of parameter space at 

which all salient constraints from neutrino physics, cosmology, and baryogenesis are satisfied 
(represented in fig. d by the points at which a given ribbon intersects the grey, vertical line 
of constant 6). At each of these points, the value of (x) is obtained by requiring that the 
neutrino masses given by equation ((J)) satisfy the experimental bounds in equation . 
For = 10^'^ GeV and 5 = me/mf^, the strong washout case (with / = 8.3 x 10~^) 
corresponds to {x) ~ 5 GeV; in the drift-and-decay case (where / = 1.5 x 10~^), (x) ~ 
50 TeV. In general, for the strong washout regime, the value of {x) increases with decreasing 
; for the drift and decay case, it decreases with decreasing . These trends continue 
until the two points converge into one at ~ 10^'' GeV, for which (x) ~ 1 TeV. 

In CHDL, the flavor structure of neutrinos rises from a hierarchy among the M<j>. in which 
6 is small. As is evident from fig. [3 it does not seem possible simultaneously to obtain a 
realistic neutrino spectrum and obtain the correct baryon number when -/Vf$^ < lO^'^ GeV. 
Since for higher we require Tr > lO^'' GeV, the reheating temperature constraints 
discussed in section |21 are of genuine concern. In particular, the constraints from BBN and 
nonthermal LSP production make things very difficult for CHDL when ^ 10^ GeV 

and BBN constraints come into play. While it is possible to get Dirac neutrinogenesis to 
work when 10^ GeV < "^/3/2 ^ 10^ GeV, this requires one of two things: either uilsp is 
light enough that the naive reheating temperature bound permits a reheating temperature 
above lO^'^ GeV (see fig. El), or else the ratio must be large enough that LSP 

annihilations are effective (fig. While it is possible to get leptogenesis to work in the 
former case, it is far from optimal in the sense that nonthermal (or non-LSP) dark matter 
is required. In the latter, thermal CDM is still possible, but the ratio m3/2/iT^LSP has to be 
be far larger than the (5/gx that results from the one-loop splitting given in equation (^. 
As discussed earlier, this means there will be great difficulty getting CHDL to work in 
loop-split SUSY when is on the PeV scale. There is, however, a window where ~ 
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Figure 7: Bands m. f — b parameter space for which the final baryon number falls within 
the range permitted by WMAP, for different choices of M$j. Here / parameterizes the 
couplings of the fields (see eq. (|54j) l. The configuration of the left (right) panel is the 
one marked with a dot in the left (right) panel of fig. El in which = 95 (03 = 4.5). The 
shaded vertical lines show the constraint on the value of 6 coming from neutrino mixings 
and masses {6 = me/m^ ~ 4.83 x 10~^ in the left and 5 = 10~^ in the right). Note that 
there are two points that yield a realistic value of Bj for a given 6 (the points at which 
the grey vertical line intersects a given ribbon). When / is small enough, the initial baryon 
number generated is just enough to be consistent with WMAP, while the washout effect 
of inverse decays and 2 <-> 2 processes is negligible. In this situation the baryon number 
generated is independent of M$ and its the final value is proportional to the CP violating 
parameter e (see eq. (Unj). The dark dotted curve shows the band of consistent baryon 
number calculated in this "drift and decay" limit. Each "X" marks the point in which 
To/H = 1 for each different M$^. At these points, the "strong washout" regime starts. 
The now active washout processes reduce an initial surfeit of baryon number (due to a larger 
/) down to an acceptable level (see fig. IBJ- 



10 GeV and misp ^ 300 GeV in which CHDL can be made to work in such scenarios. For 
7713/2 ^ 10^ GeV (where the gravitino decays rapidly and the LSP relic abundance is entirely 
thermal), reheating temperatures greater than lO^'' GeV are permitted, and CHDL succeeds 
in providing an explanation for the origin of the observed baryon asymmetry and neutrino 
mixings. Since constraints from gluino cosmology become relevant for ^ 10^'' GeV |33j . 
it appears that split supersymmetry models where 10^ GeV < ^ ^^^^ GeV are the 

"natural habitat" , so to speak, for thermal Dirac neutrinogenesis with thermal CDM. 
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Figure 8: These two plots show the evolution of baryon number B for = 10^^ GeV 
and b = nie/m^ = 4.83 x 10~^ (the CHDL value), in the two different regimes that produce 
a realistic value for the final baryon number of the universe, Bp. For / = V A23/123 = 
1.5 X 10"'^, as shown in the left panel, the effects of 2 <-> 2 lepton-number-changing processes 
are negligible and the final baryon-to-entropy ratio is the same as that initially produced by 
(j) and (p decays. For stronger coupling / = 3.8 x 10^^, as shown in the right panel, baryon 
number is initially overproduced, but 2 2 processes, which are stronger for stronger 
coupling, reduce B to an acceptable level by the time they freeze out. The left and right 
panels correspond respectively to the lower and upper parts of the 'ribbon' in fig. [3 



6 Conclusion 



Dirac neutrinogenesis is an interesting alternative to the standard Major ana lepto- 
genesis, and one which is equally capable of explaining both the baryon asymmetry of 
the universe and the smallness of neutrino masses (without the seesaw mechanism). We 
have shown that thermal Dirac neutrinogenesis is not only an interesting theory, but a 
genuinely viable phenomenological model: it is simultaneously capable of producing a 
baryon-to-photon ratio r] for the universe that matches that observed by WMAP, yielding 
a neutrino spectrum that agrees with the predictions of solar and atmospheric neutrino 
data, and avoiding reheating temperature constraints and other potential cosmological and 
astrophysical problems. Furthermore, it is capable of doing all this in the context of a 
simple setup in which the observed neutrino masses and mixings are explained by imposing 
(anti) symmetry conditions on couplings in the superpotential and a hierarchy among the 
masses M$- of the additional heavy fields in the theory. 

We have also shown that Dirac neutrinogenesis is naturally accommodated in split 
supersymmetry regimes involving either a heavy gravitino (with a mass above around 
10^ GeV) or a very large hierarchy between the gravitino mass and the mass of the LSP. 
In models with a lighter gravitino and only a small splitting between 771,3/2 and misPy LSP 
annihilations are generally ineffective, the reheat temperature Tr associated with inflation 
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cannot be raised above ^ 10 GeV, and constraints on become quite severe. As 

a result, leptogenesis becomes extremely difficult (for ease of reference, we classify the 
different gravitino-mass regimes with regard to leptogenesis considerations in table 12) • It 
is interesting to note that Dirac neutrinogenesis, at least in its constrained hierarchical 
form (CHDL), can be realized in simple, one- loop AMSB in a small region of parameter 
space in which M<j>j ~ 10^'' GeV and rriLSP ^ 300 GeV (which sets 1713/2 ~ 10^ GeV). In 
more general split SUSY models, cosmological constraints are less severe. In models where 
^ 10^ GeV, the constraints altogether disappear and Dirac neutrinogenesis can be 
made to work without sacrificing the possibility for thermal LSP cold dark matter. Thus 
Dirac neutrinogenesis, at least in a regime with heavy gravitinos, provides an interesting 
and viable alternative to Major ana leptogenesis. 



Gravitino Mass 
Range (GeV) 


Maximum 
M$, (GeV) 


Workability 
of CHDL? 


Comments 


"13/2 < 10"^ 




Very Low 


G decay during or after BBN. 
Insufficient 77 generated. 


10^ < 1713/2 < 10* 


10° - lOlO 
or higher 


Depends on 
tni^yp and the 
ratio 


LSP annihilations ineffective unless 
"^LSP 1 '^'i jl small. constrained 
by nonthermal LSP abundance from G 
decay. Loop-split SUSY works only for 
"^3/2 ^ GeV. More general split 
SUSY theories can be successful. 


10* S "13/2 £ 10^° 


None 


Excellent 


^l^SP thermal, since G decays 
before LSP freeze-out. 

A/^^ > 10"'""'" GeV allowed, v spectrum 
requirements compatible with CHDL. 


10^° < m3/2 < 10^^ 


None 


Questionable 
(depends on 

gluino 
properties) 


V sector and baryogenesis okay, but 
model may have a cosmological gluino 
problem. 


10^^ < "13/2 


None 


Very Low 


Potential gluino problem becomes a 
serious concern. 



Table 2: The various gravitino mass regimes for split supersymmetry models and the 
viability of thermal Constrained Hierarchical Dirac Leptogenesis (CHDL) in each case. 

Another asset of the theory is that there are experimental checks on its viability. The 
major prediction of Dirac neutrinogenesis is that neutrinoless double-beta decay will not 
be observed to any degree, for this process relies on the existence of a Majorana mass term 
for right-handed neutrinos. The discovery of such a process experimentally would rule the 
theory out. In addition, forthcoming results from MiniBooNE should either confirm or deny 
the LSND result, which will reveal whether or not the neutrino spectrum produced by Dirac 
neutrinogenesis is in fact the one present in nature. Finally, thermal Dirac neutrinogenesis 
yields some definite predictions about the neutrino mixing parameter sin^ia (see fig. El, the 
value of which will be measured in future experiments. 
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A Derivation of the Boltzmann Equations 



In this appendix, we derive the Boltzmann Equations that appear in equations (|63j) - 
H62|) . following the methods of ;46j. The Boltzmann equations for any particle species a in 
the early universe can be written in terms of 1^ = Ua/s, where Ua is the number density of 
a and s is the entropy density, as j24j 



dYA 
dt 



1 f d?Pa d?Pi d?Pj d?Pk 



s J (2^)3 (2^)3 (27r)3'"(27r)3 



n=a,i,j,...k 

^ \M{a...i ^ j...k)\\fa...fi) - \M{j...k ^ a...i)\\f,...fk) 



int. 



int. 



(71) 



where fi is the phase-space distribution function of particle i, \Ai{a...i — > j...k)\'^ are the 
squared matrix elements for particle-number-changing interactions involving a, and the 
sums are over all interaction processes which create or destroy a. Since for the scalar fields 
(j) and (p'^, the leading processes include only decays and inverse decays, the Boltzmann 
equations for the abundances and Yr are given by 



dY^ 



1 . 



dt 



-A 



12 



(72) 



and 



dt 



-A 



12 



(73) 



where A|' '^ has been used as a shorthand to denote the appropriate phase-space integral 
and we use the notation / and to denote the usual Lepton and Higgs-up doublets of 
SU{2). 

Let us also define the particle asymmetry A = (Ya — Y^), where Ya and Y^ are the 
entropy normalized abundances of particle a and its conjugate a^. For a relativistic species 
a (i.e. T ^ nia) and as long as both A and /i^ are small, we have A ~ ^ga{e^^°'^'^ — 1)) 
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1 + :;r-^ where Qa is the number of degrees of freedom of particle a, 



which leads to e'^"/'^ ~ 1 + 
s is the entropy density, and is the photon density. If we apply the energy conservation 
relation to the inverse decay processes of the form (a + i — > 6) appearing in equations (|72|) 
and H73|l . in which particle i is in chemical equilibrium with other particles in the thermal 
bath (i.e. /ij = 0), the previous relation allows us to write 



(74) 



where /a'' is the equilibrium distribution of a, and the negative sign occurs when a is a 
conjugate particle. 

It will also be useful to relate the particle asymmetries of Leptons to the net Lepton 
number abundance carried by each Lepton species. We define the Lepton abundances as 













= Y~-Y~ 







Relation (|74() . along with the CPT conservation relation 

\M{a^ij)\' = \M{i''f 

allows us to simplify the Boltzmann equations for (p and 
be stated in terms of Li, L^,^, L-p^ and as 



,cm2 



(75) 
(76) 
(77) 
(78) 



(79) 



significantly. The results may 



dY^ 



dt 

dYA,c 



dt 



(27r)3 
(27r)3 



eq 



eq 



where the interaction rates Tl 
the relations 



t> ^ i + x) and = r 

d^Pi d'^P' 



(2^)3 (27r)3 
(fpi d^pi 
(2^)3 (27r)3 



\M{(t> 
|A^(0 



(80) 
(81) 

z^^ + H^) are defined by 

(82) 
(83) 



with r^, R being the rates of the conjugate processes and F^i = L^ + Tr is the total decay 
rate for (/>, etc. given in ((T^ . Because of supersymmetry, the total decay rate of the 
fermion components of the heavy supermultiplets <I> and $ will also be F/), with the same 
and Tr as their partial rates. 



We also define the Lepton asymmetry Lw 



Y4, — Y^c associated to the new heavy 



leptons. After integrating over the incoming momentum 

dLs, 



we obtain for La 



dz 



-(Fd) L 



s 
2n 



-Yl'e 
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eq 



Lj 



{Tl) + L,^{Tr)] (84) 
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Here, (Ta) is the respective decay rate averaged over time-dilation factors 



K,iM^/T) 
^^^^ = K,{M^/T) 

where Ki{x) and K2{x) are modified Bessel functions. The Boltzmann equations for the 
scalar component of the $ superfield (and its conjugate), as well as the ones for the fermion 
components of $ and are obtained in a similar manner. 

We now turn to address the evolution of the Lepton number abundance of the particle 
species vr, i, and in which we must take into account the effect of 2 ^ 2 processes 
which transfer lepton number between Li, Ly^, Lj, and L^^^. We will make one important 
distinction among them: the rates for interactions which shuffle lepton number between 
£, £, vr, and the right-handed charged lepton and sleptons fields cr and cr will be much 
larger than those for the interactions which shuffle lepton number between ur and any 
of these other fields. This is due to the fact that vr only interacts via the processes 
pictured in figure which involve a virtual (f)^, 4>—, etc. while all of the other fields 
either take part in SU{2) and/or U{1)y gauge interactions, or — in the case of i/r — left- 
right slepton equilibration through an assumed large {F^) term. We will represent the 
effects of these rapid equilibration processes by including terms Tia (where A is the relevant 
particle asymmetry) to represent them in the Boltzmann equations. The slower 2 <-> 2 
processes through which right handed neutrinos vr interact with lepton doublets i and 
i and with right handed sneutrino Vr will be included separately as C^^^^, C^^^j and 
Cy^^y^. To simplify further the notation we will define the terms Fa to account for the 
collective contribution from decays (and inverse decays) of the fermionic components of $ 
and <I> (which we will not write explicitly, being of similar form to the contribution from 
the scalar components). 

We obtain for the Lepton number abundance L^: 

^ = -^{'^D){Y,.+Yll)+L,{T,)-Yll^L,({T:^)-]-e{T 



+Fe + ^e + C^^^^ (86) 
Similarly, the equations for Ly^^, L^, and Lp^ are given by 

'-^ = 'i^o){n,+Y;i)-L,jT„) + Y;i±L.,{ir,)-'-.r,) 

+ Fyj^ - C^^^J - Cyj^^yj^ " C^^^^ (87) 



^ = ^(rz.)(y,^ + y;!)-L,^(r.)-y-^L,4(r,) + i.r. 

+ Fyj^ + ^Vr + Cy^^y^. (89) 

Before going any further we still need to compute the terms C^^^j, Cy^^y^ and C^^^^ 
corresponding to the 2^2 processes mediated by heavy fields. 
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We begin by calculating C^^^^, which is given by 

[\m'{h:i.% Ix)? - \M'{lx - ^'kK)?) 



C ~ = AA34e-(i^i+i^2)/T 



here refers to the amplitude for the specified 2^2 process to which we have substracted 
the contribution from the resonant intermediate state (RIS) in which a real field (j)^ is 
produced and then decayed into the 2 particle final state. The RIS contribution must be 
substracted since we have already counted contributions from decays of real (j)^ fields. 

The leading term in the difference between a 2 <-> 2 process of the form ab — > ij involving 
a heavy intermediary k and its conjugate process depends on the contribution of the on-shell 
(resonant) intermediate state: 

\M\ah ^ ij)l' - \M\ij ^ a6)|2 = \Mnxs{ij ^ ab)\'' - \Mnis{ah ^ (91) 

with 

\Mnxs{ah ^ ij^ ^ -^Ks - nil)\M{ah ^ k)\^ x \M{k ^ (92) 

where k represents the intermediate-state particle, and s is the usual kinematic variable 
s = (pI"- +^2")^- In our case, making use of the equality T^r^j — T/^r^ = we find 

\Mnxs(ix - ^'rK)? - \Mnxs{K^'k - ~^x)? ^ e^^5{s - ml)\Mit\' (93) 

and substituting this result in eq. (|9U|) we finally obtain 

C^^^j = 2eY;i{Tn) + {L.,-\l^) n,{va^^^j+vaj^J. (94) 

For the other sets of 2 •s-^ 2 processes, C^^^i,^ and C^^^^i, the procedure is essentially the 
same. The rates (T^^^^) = n^{va^^^^ + va^^^J, {T^j^^e) = n^ivauj^^e + vag,^^^) and 
(Ty^^y^) = n-y{vayj^^'Pj^ + va-c^^y^) , associated with these interactions are calculated in 
section |SJ We will denote the total contribution from these processes as (r2^2)- 

At this point, the full Boltzmann system comprises fourteen individual differential 
equations: four for (p^, (/>_, and their conjugates (or alternatively, the abundances ^(/i^, 

, etc.); four for the the Lepton asymmetries L^, Ly^, Lj, and Lp^; two for the Lepton 
asymmetries = — e^j and Lg^ = — of the right-handed charged lepton 
and slepton fields; and an additional four for the fermionic superpartners in the $ and 
$ supermultiplets. 

Let us make the assumption that the rates for the processes denoted by 'Sa are suffi- 
ciently rapid that chemical equilibrium is achieved among the particle species that interact 
via these processes. In this case, any lepton number stored in i, £, i>R, sr, or e/j should 
be rapidly distributed among all of these particles in proportion to the relative number of 
degrees of freedom for each field. The lepton number stored in each field then becomes 

Li,Lj Leg, (95) 

, Lej^, L^jl (96) 
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where we have defined an equihbrium lepton number L^q which represents the effects of rapid 
equilibration through Ha- Because of this rapid equilibration, it is no longer necessary to 
keep track of the individual lepton numbers L^, L^, Lj;^, Lg^, and Lg^: we need only one 
differential equation representing the total lepton number in this equilibrated, aggregate 
sector Lagg = + L^^ + Lj+ Lg^ + Lg^ = ^L^q- The equation is 

^"^^^^ - <^D)iYl-YlA+{T,)(u+LA+{T^)U 



dt ' ^' \ 



~^Zi;^^3^{^^^^ + ^^^0 + " 7^"^^^ (r2^2) +Fnet, (97) 

where Fnet represents the overall contribution to Lagg from the fermionic components of $ 
and and small terms proportional to e times Li, L^^, ^i^^ ^ur have been dropped. Here, 
we have used the fact that the Boltzmann equations for Lg^ and Lg^ are trivial, consisting 
of only terms, which all cancel after taking the sum of all Lepton abundances. 

The one task that remains is to deal with the fermionic components of <I> and <1>. However, 
it can be shown that in the approximation that the T,a interactions are rapid, the decay and 
inverse decay contributions to Lagg and L^^^ from these fermions are the same as those from 
the scalars, and thus no additional Boltzmann equations for ■0$, etc. are needed. We 
therefore require only five differential equations to describe the evolution of lepton number 
in the universe. These are 

^ = -2e{To){Yl^-Y;i) + {r,)iL,^+L,^) + {rn)L,^ 

-2Lagg ((Td),^ + (Tl),^) + (L,^ - ^Lagg){r2^2) (98) 

^ = 2e{Tn){Yl-Y;i) + L,jT^)-2L,jTn)^, 

-{L,^-^Lagg){r2^2) (99) 

dY^ 1 1 

= -{Tom^ - Y;1) + -LaggiVL),, + -L,, (r^),, (100) 

dL^ 

= -{rD)L.+2Lagg{rL),^ + 2L,^{rR),^ (loi) 



dt 
dLj, 



{rD)L4,_ + 2Lagg{rD)jn, (102) 



dt ' ^' 

eq eq eq 

1 1 ^0 

where {Td)io = n — -{^d), = 7 — -{^l) and (Tr)^^ = — -{Td)- Once again, 

small terms proportional to cLagg or eL^^ have been neglected. Of course we have not yet 
included the coupling of baryon number B to Lagg, which is accomplished by introducing 
a differential equation for B (making a total of six coupled equations) and the addition 
of a sphaleron interaction term which mixes B and Lagg. With this addition one obtains 
equations ^ - (|H7|) . 

Before closing this appendix we can easily solve these equations in the "drift and decay" 
limit, in which we assume that all the rates {Ta) are much smaller than the rate of expansion 
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of the universe H. In that hmit the equations simphfy greatly: 

dL 



from which we can write 



^ = MrD){Yl-Y;i) (104) 



AT AT dV^ 



2e — ^. (106) 

dt dt dt ^ ^ 

Since the inital values of L^^ and Lagg vanish, the final abundances after a long time will 
be 

Lir' = -L{:r' = 2<! (t = o) (lo?) 



where we have assumed that yj^"'"'^ = and y^"***«' = y^"? (t = 0) with t = defined as the 
time in which T = M$. 
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